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ABSTRACT 


A weakest  link  theory  for  failure  of  brittle  materials  containing  randomly 
oriented  cracks  uniformly  distributed  throughout  the  volume  is  compared  with 
an  analogous  theory  for  surface  distributed  cracks.  In  the  latter  theory  all 
crack  planes  are  assumed  to  be  normal  to  the  surface  in  which  they  are  located. 
In  both  theories  failure  is  assumed  to  occur  when  the  component  of  stress 
normal  to  a crack  plane  exceeds  the  strength  of  the  crack. 

It  is  shown  for  both  theories  that  the  distribution  of  crack  critical 
strength  can  be  given  in  closed  form  when  the  failure  statistics  in  simple 
tension  are  known.  It  is  also  shown  that  the  predicted  biaxial  failure  sta- 
tistics based  on  given  uniaxial  statistics  are  the  same,  regardless  of  which 
theory  is  used.  Evidence  is  cited  which  suggests  that  the  latter  result  may 
hold  for  shear-sensitive  as  well  as  shear-insensitive  cracks. 


i 


9 


' 

• • > 
ivEfl  .... 


» 


t 


t 


TABLE  OF  CONTENTS 

PAGE 

1 

I .  INTRODUCTION  

II.  EVALUATION  OF  CRACK  DENSITY  FUNCTION  

III.  EQUIVALENCE  OF  THE  TWO  THEORIES  

12 

IV.  DISCUSSION  

14 

APPENDIX  

17 

references  


f 


» 


ii 


...  .-** ' 

;/*  • " ■’> 


I . INTRODUCTION 


Many  structural  materials,  especially  those  with  high  strength  at  high 
temperature  and  those  that  are  transparent  to  electromagnetic  radiation,  tend 
to  be  brittle.  Brittle  materials  exhibit  a substantial  dispersion  in  fracture 
stress  which  must  be  taken  into  account  in  design.  This  is  done  by  making  use 
of  the  measured  or  theoretically  determined  fracture  statistics  for  the  material 
under  the  stress  conditions  encountered  in  service. 

The  statistical  fracture  theory  most  frequently  employed  is  that  proposed 
by  Weibull ^ in  1939.  This  theory  assumed  that  the  variation  in  measured  strength 
of  nominally  identical  specimens  is  due  to  the  presence  of  invisible  flaws. 

The  flaws  have  a distribution  in  strength  characteristic  of  the  material.  As 
a result,  different  volume  elements  have  different  flaw  contents  and  different 
strengths.  A structure  is  assumed  to  fail  when  the  weakest  element  is  over- 
loaded, just  as  a chain  breaks  when  the  strength  of  its  weakest  link  is 
exceeded.  Weibull  worked  out  the  fracture  statistics  of  the  weakest  link 
concept,  and  in  addition  gave  without  proof  a procedure  for  finding  the  frac- 
ture statistics  for  polyaxial  stress  states  when  they  are  known  for  uniaxial 
tension.  Weibull* s rules  are  very  simple  to  apply  when  the  uniaxial  data  are 
represented  in  terms  of  Weibull *s  2-parameter  formula. 

2 

A new  weakest  link  theory  was  proposed  recently  by  Batdorf  and  Crose. 

This  theory  differed  from  Weibull *s  in  several  respects:  (1)  it  identified 

the  flaws  as  flat  cracks,  an  assumption  which  implies  that  flaws  have  a 
strength  which  varies  with  their  orientation  with  respect  to  the  applied 
stresses,  (2)  it  generalized  Weibull's  3-parameter  representation  of  experi- 
mental data  and  made  more  accurate  fitting  possible  by  the  use  of  a Taylor 


expansion  with  an  arbitrary  number  of  coefficients,  and  (3)  it  assumed  that 
only  the  component  of  stress  normal  to  the  crack  plane  contributes  to  the 
fracture  of  the  crack.  The  latter  assumption  is  a fair  approximation  and  has 
the  advantage  that  use  of  it  avoids  the  necessity  of  specifying  crack  shapes 
and  Poisson's  ratio.  The  theory  was  formulated  for  volume  distributed  cracks 
in  a macroscopically  isotropic  material,  so  all  crack  orientations  were  con- 
sidered equally  likely. 

An  analogous  theory  for  surface-distributed  cracks  was  subsequently 
3 

developed  by  Batdorf.  This  theory  assumed  that  the  crack  plane  is  always 
normal  to  the  surface.  Thus  the  orientation  of  the  surface  can  be  specified 
with  a single  parameter,  whereas  the  orientation  of  an  interior  crack  requires 
two  parameters. 

In  References  2 and  3 the  coefficients  in  the  Taylor  series  used  to  repre- 
sent the  distribution  of  crack  strengths  were  determined  by  inverting  a matrix. 
This  implies  restricting  the  series  expansion  to  a limited  number  of  terms, 
and  thus  involves  an  approximation.  It  was  subsequently  shown  by  D.  J.  Chang 
that  an  integral  equation  approach  can  be  used  which  is  in  principle  exact. 

In  the  following  section  exact  approaches  are  given  for  both  volume  distributed 
and  surface  distributed  cracks.  Also  a rather  surprising  fact  relating  the  two 
theories  is  proved.  This  is,  if  the  statistics  of  biaxial  tension  are  derived 
from  known  statistics  of  uniaxial  tension,  the  results  obtained  are  the  same 
irrespective  of  which  theory  is  used.  As  a result  the  simpler  surface  crack 
theory  can  be  used  to  find  the  fracture  statistics  of  biaxially  stressed 
structures  with  volume  distributed  cracks . 


II.  EVALUATION  OF  CRACK  DENSITY  FUNCTION 
Volume  Distributed  Cracks 

According  to  the  statistical  fracture  theory  outlined  in  Reference  2, 
the  probability  of  survival  of  a volume  V of  material  subjected  to  a uniform 
stress  state  L is 


P (L)  = exp 
s 


1 SUL, a ) dN(a  ) 
cr cr 

4tt  da 

cr 


(1) 


In  this  expression  a^  is  the  maximum  principal  tensile  stress.  N(ac^_)  is  t*le 

number  of  cracks  per  unit  volume  with  a critical  stress  less  than  or  equal  to 

a , and  ft  is  a solid  angle  such  that  a crack  with  critical  stress  O will 
cr  cr 

be  fractured  if  and  only  if  its  normal  lies  within  the  solid  angle,  ft  can  be 

calculated  for  any  stress  state,  so  the  corresponding  probability  of  failure 

can  be  determined  if  N(a  ) is  known.  N(a  ) can  be  evaluated  if  P (a)  is 

cr  cr  s 

known  for  any  stress  state.  Generally  the  available  data  are  for  uniaxial 

tension.  In  Reference  2 an  approximate  technique  was  employed  for  obtaining 

N(a  ) from  data  in  simple  tension,  which  involved  inverting  the  matrix  of  an 
cr 

ill  conditioned  set  of  simultaneous  linear  algebraic  equations.  A better  tech- 
nique is  to  treat  Equation  (1)  as  an  integral  equation.  An  elementary  method 

4 

of  doing  this  suggested  by  B . Budiansky  follows . 

It  is  shown  in  Reference  2 that  for  the  case  of  simple  tension,  the 
assumption  that  only  the  component  of  stress  normal  to  the  crack  plane  con- 
tributes to  fracture  leads  to 


4ir 


1 


(2) 


3 


Thus, 


p (a)  = exp 
s 


ro  / /o- ' 

K (W- 


a \ dN(a  ) 
— ' cr 


do 


do 


cr 


cr 


Integrating  by  parts  and  making  use  of  the  fact  that  N(0) 
term  vanishes  and 


(3) 


= 0,  the  integrated 


P = exp 
s 


/•O  N(0  ) 

_ v / do 

J0  2^  Cr 

cr 


(4) 


Thus, 


>/°toPs“"I  L 


o N(o  ) 


cr 


do 


0 \ a 


cr 


h *n  ps) 


V N(o) 
s/  2 


(5) 


(6) 


from  which 


VN(O)  = ~ 2 jA  £n  Psj 


(7) 


Surface  Distributed  Cracks 

According  to  the  analogous  theory  for  surface  cracks  of  Reference  3,  the 
probability  of  survival  in  any  stress  state  Z is  given  by 


P (E)  = exp 
s 


dN(o  ) 

_ A / “ 2L-  do 

J0  77  d0cr  cr 


(8) 


where  A is  the  area  and  u)  is  an  angle  such  that  a crack  of  critical  stress  Ocr 
will  be  fractured  if  its  normal  lies  in  to.  It  is  shown  in  Reference  3 that  for 

uniaxial  tension 


P (a)  = exp 
s 


— f 

" 17  Jo 


dN  -1,  / cr  Jn 

da  Vo  cr 

cr 


(9) 
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Taking  the  logarithm  of  both  sides  and  integrating  by  parts  we  obtain 
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III.  EQUIVALENCE  OE  THE  TWO  THEORIES 


In  this  section  we  assume  that  the  probability  of  survival  is  known  for 
simple  tension  and  calculate  the  probability  of  survival  for  biaxial  tension. 
The  latter  result  turns  out  to  be  the  same  for  surface  distributed  cracks  as 
it  is  for  volume  distributed  cracks. 

If  the  probability  of  survival  in  simple  tension  is  a reasonably  well- 
behaved  function  of  the  stress  level,  it  can  be  expressed  to  arbitrary  accuracy 
by  the  power  series 

Pg(a)  = exp  [-  E a±  cf1]  . (13) 

Using  the  theory  of  the  previous  section,  we  use  equation  (13)  to  find 
VN(CTci_)  and  AN(°cr) > respectively,  and  then  calculate  the  probability  of 
failure  for  a biaxial  stress  state. 

Volume  Distributed  Cracks 


Substituting  Equation  (13)  into  Equation  (7)  we  obtain 


VN(a) 


2 !^a r ai °i  • 


Carrying  out  the  indicated  operations  we  obtain 


VN(o  ) = 2 E (i  + 0. 5)a  a*  . (15) 

cr  ^ i cr 

For  each  of  the  special  cases  of  simple  tension  and  equibiaxial  tension, 

S7/4tt  is  a simple  function  of  a /a.  For  other  cases  we  can  define  fi  as 
r cr 


■r< 


sin  9 d 0 , 


when  the  coordinate  system  used  is  shown  in  Figure  1.  The  limits  <p  , 0^ 
and  ©2  are  chosen  in  such  a manner  that  the  integral  is  taken  only  over  the 


- . J • 

<1'»  *>  I .<K 


range  of  angles  and  0 within  which  0^  > a^.  Alternatively,  we  can  define 
f 1 by 

2tt  /.tt 

d p I sin  9 H(a  , a ) d8  , (17) 

Jo  cr  n 

where 

H(a  , a ) = 1 for  a < a (18a) 

cr  n cr  n 

= 0 for  a >o  . (18b) 

cr  n 

Substituting  Equation  (17)  into  Equation  (1)  and  integrating  over  Ocr  we  obtain 

P (Z)  = exp 
s 


- V 


/•2tt  /.tt 

J d*J 

J0  J0 


d0  sin  0 


N(a  ) 
n 


(19) 


We  note  in  passing  that  this  equation  is  equivalent  to  the  procedure  for  poly- 
axial  stress  states  given  without  proof  by  Weibull  in  Reference  1. 

Substituting  Equation  (15)  into  Equation  (19)  and  making  use  of  the 
relation 


2 2 2 2 2 2 
o = (a, cos  P + a„  sin  p)  sin  0 = a,  (cos  P + Ksin  <p)sin  0 
n 1 2 1 

where  K = a2^ai  * 

we  obtain 

1 i 2 2 1 

to  P =~  118.(1  + 0,5)o.  I (cos  P + K sin  P)  do 

3 277^^1  1 •'o 

x f sin2l+1  0 d 0 
J0 

= - — £ ai(i  + 0.5)0^  4>(i,  K)  B(i  + 1,  y)  , 


(20) 

(21) 


(22) 


where  B is  the  beta  function  and 


The  substitution 


sin  'P  - z (24) 

leads  to 


If  we  let 


then  Equation  (27)  becomes 


. rl  i-1 

= Z i a a 1 / ?~= 

i 1 cr  *A)  y i - 


AN(a  ) = 
cr 


i r(1)rlf 

= E i a a 7— ^ 

i 1 cr  r/i+  \ 


= E a.  a , 

i cr  * 


where 


ai  = ai 


ra+i)  r(j) 

1 rK) 


In  the  case  of  biaxial  tension.  Reference  3 gives  for  the  probability  of 


survival 


..a/1 

IT 


20  da 

_ cr  da  cr 
0 cr 


where 


. / a -a 

0 „ = cos-1  \ 1 - — /y  ~ for  Ka  < a < a. 
cr  \J  (1— K)  1 cr  1 


(34a) 


for  a < Ka, 
cr  1 


(34b) 


Integrating  by  parts 


TT  J„ 


N(a  ) da 

cr  cr 


Ka,  / a,  - a / a,  - a 

V V1  - V1  - K> 


Making  the  substitution 


0,-0 

1 cr 

ct1(1  - K) 


and  substituting  Equation  (31)  into  Equation  (35)  we  obtain 


, yl  £a  a dz 

tnP  --\(  -yr-= 

s ’•'o  yr  yr^ 


Now  from  (36)  it  follows  that 


a = o,  - z(l-K)o1 
cr  1 i 


Substituting  Equations  (38)  and  (32)  into  Equation  (37)  we  obtain 


r(i+i)rlf  J 

Jin  P = a,  (a.)1  -7- »(i,  K)  , 

s 27T  1 1 T|  i'+  jj 


which  is  identical  to  Equation  (26)  . 


IV.  DISCUSSION 


What  we  have  just  demonstrated  is  that  when  the  uniaxial  stress  statistics 
are  known,  the  biaxial  stress  statistics  turn  out  to  be  the  same,  irrespective 
of  whether  we  employ  the  theory  for  volume  distributed  cracks  or  the  theory 
for  surface  distributed  cracks.  Since  the  statistics  for  any  stress  ratio  are 
uniquely  related  to  the  statistics  for  any  other  stress  ratio,  it  is  not  neces- 
sary to  use  uniaxial  stress  statistics  as  the  starting  point.  The  statistics 
for  any  convenient  stress  ratio  can  be  taken  as  known  and  the  statistics  for 
any  other  stress  ratio  can  be  calculated;  using  either  theory,  the  same  result 
will  be  obtained.  Since  surface  theory  involves  an  integration  over  only  one 
variable  whereas  volume  theory  requires  integration  over  two,  it  is  generally 
advantageous  to  use  surface  theory. 

The  preceding  remarks  apply  to  the  case  of  shear-insensitive  cracks,  i.e., 
to  the  case  in  which  only  the  normal  stress  contributes  to  fracture.  It  is 
known  that  the  shear  acting  on  the  crack  plane  also  contributes  to  fracture, 
and  it  would  therefore  be  of  considerable  interest  to  know  whether  the  same 
result  holds  for  shear-sensitive  cracks.  This  is  not  an  easy  question  to 
settle  in  complete  generality.  First,  there  is  as  yet  no  concensus  in  the 
fracture  mechanics  community  as  to  how  to  treat  mixed  mode  fracture.  Thus,  we 
would  have  not  one  but  several  theories  of  shear-sensitive  cracks  to  investigate. 
Second,  even  the  simple  case  of  shear-insensitive  cracks  required  a somewhat 
involved  theoretical  analysis,  and  it  is  to  be  expected  that  the  analysis  for 
shear -sensitive  cracks,  in  some  cases  at  least,  would  be  considerably  more 
involved.  However,  it  is  easy  to  show  that  the  same  rule  applies  in  a partic- 
ularly simple  special  case.  This  is  the  case  in  which  Pf(0,  a)  is  calculated 
from  P^(a,  0)  when  the  effective  stress  leading  to  fracture  is  assumed  to  be 


i 


where  cr  and  x are  the  normal  and  shear  stress  acting  on  the  crack  plane 
respectively.  A proof  for  this  statement  is  outlined  in  the  Appendix. 

Although  the  matter  is  not  yet  settled,  it  is  the  present  view  of  the  writers 
that,  at  least  for  some  shear-sensitive  crack  models,  the  same  rule  applies. 

Finally,  it  should  be  noted  that  the  equivalence  between  volume  and 
surface  distributed  crack  theories  is  restricted  to  the  question  of  using  the 
known  fracture  statistics  for  one  stress  ratio  to  calculate  the  statistics  for 
another  stress  ratio.  The  distributions  of  crack  critical  stress  deduced  from 
a given  set  of  data,  i.e.,  the  functional  forms  found  for  11(0^),  will  be  quite 
different  for  the  two  theories. 


APPENDIX 


Let  us  consider  a particular  type  of  shear -sensitive  crack,  i.e.  that  for 
which  fracture  occurs  when  the  function 


reaches  a critical  value.  Since 

2 2 

a = a,  cos  0 + sin  0 
n 1 z 

T = (cr^  - a2)sin  0 cos  0 , 


(Al) 


(A2a) 

(A2b) 


it  follows  that  for  uniaxial  tension 

0eff  = 0 COS  6 * (A3) 

For  the  volume  distributed  crack  theory.  Equation  (1)  still  applies,  but  ft 
is  changed: 


ft_ 

4tt 


1 


a / a 
cr 


uniaxial  tension 


(A4a) 


= \/l  - <aCIVa> 


equibiaxial  tension 


(A4b) 


When  Pg(cr,  0)  is  known  we  can  find  N(a)  in  the  following  manner.  Com- 


bining Equations  (1)  and  (A4a), 


£n  P = - V f (i  - N'  (a  )d a 

s ' a / cr  cr 


(A5) 


Integrating  by  parts  we  obtain 

a 

0 

whence 


£n  P 


■ if 

0 Jn 


N(o  )da  , 
cr  cr  ’ 


CA6) 


v N(o)  = - [o  £n  P Ca,0)l 


(A7) 


The  plan  of  attack  here  is  to  assume  a very  general  form  for  P (a,  0)  and 

8 

use  volume  distributed  crack  theory  to  obtain  P8(a,  0) . Using  this  value  of 

Pg(o,  0)  as  given,  we  will  use  surface  distributed  crack  theory  to  calculate 

P (0,  0),  and  show  that  it  is  the  same  as  what  we  used  at  the  outset, 
s 

With  this  in  mind  we  assume 

In  P s(0,  0)  = -E  a^a1  . (A8) 

Combining  Equations  (1),  (A4b) , (A7) , and  (A8) , we  obtain 


In  P (a,  a)  =-E  a (i  + l)i  f 
8 i 1 J0 


f o1-1 

Jn  cr 


1 - (o  /a)  da 

cr  cr 


Substitution  of 


x = (acr/o)‘ 


(A10) 


converts  the  integral  into  a standard  form  for  the  beta  function,  and  the 
result  expressed  in  terms  of  gamma  functions  is 


in  P (a,  a)  = - E a. a 
s 1 


(All) 


where 


a = a (i  + l)i  I(0-5i.)£(1.5) 
i iU  l)i  r(0.51  + 1.5) 


(A12) 


Now,  according  to  the  surface  distributed  crack  theory, 


N(0)  =-£n  P (o,  a) 
s 


(A13) 


For  the  shear-sensitive  cracks  under  consideration, 


£n  P (a,  0)  =— A /N'(o  ) cos  (a  /a) do 

s J cr  cr  cr 


(A14) 


Integrating  by  parts,  and  using  Equations  (AID) , (All),  and  (A13) , we  obtain 


£n  Pg(0,  0)  = 


1 0 . 5 ( i - 1) 


i i vu,-’vl 

i E a a1  / * 

* 1 Jo  yrr 


(A15) 


i - i i r(o.5i  + o.5)r(o.5) 
TT  L i + 1 ° r(0.5i  + 1) 


(A16) 


■ <^  . •>  ■ l -M 


Substitution  of  Equation  (A12)  into  Equation  (A16)  and  some  manipulation 
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